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ABSTRACT 

The relic cosmic background neutrinos accompanying the cosmic microwave background (CMB) photons 
may hide a universal lepton asymmetry orders of magnitude larger than the universal baryon asymmetry. At 
present, the only direct way to probe such an asymmetry is through its effect on the abundances of the light 
elements produced during primordial nucleosynthesis. The relic light element abundances also depend on the 
baryon asymmetry, parameterized by the baryon density parameter (j/b = "b / n^), and on the early-universe 
expansion rate, parameterized by the expansion rate factor {S = H'/H) or, equivalently by the effective number 
of neuti-inos (N,, = 3 + 43(5^- l)/7). We use data from the CMB (and Large Scale Structure: LSS) along 
with the observationally-inferred relic abundances of deuterium and helium-4 to provide new bounds on the 
universal lepton asymmetry, finding for 77L, the analog of ?7b, 0.072±0.053 if it is assumed that = 3 and, 
0.115±0.095 along with = 3.3!*J:^, if is free to vary. 

Subject headings: lepton asymmetry — neutrinos — early universe 
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1. INTRODUCTION 

The standard models of particle physics and cosmology assume that in the early, radiation-dominated universe only the known, 
massless or light {mc^ <C kT) particles, including three flavors of light, active neutrinos (N,y = 3), contribute to energy density 
driving the universal expansion {H^ = ^TrGp/3; p = pr). It is also generally assumed that any universal lepton asymmetry is very 
small, of order the baryon asymmetry.-' In analogy with the measure of the baryon asymmetry provided by ?7b, an asymmetry 
between neutrinos and antineutrinos of flavor a {a = e, ji, r) can be described in terms of the neutrino chemical potential /Xq, or, 
in terms of the dimensionless degeneracy parameter = IJ-a/kT by. 
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For Big Bang Nucleosynthesis (BBN) the electron neutrinos play a key role through their charged-current weak interactions, 
which help to regulate the neutron-to-proton ratio {p + e~ ^ n + Ve, n + e^ ^ p + Ve, n ^ p + e~ + Ve). Since the BBN-predicted 
"^He abundance is, to a very good approximation, determined by the neutron-to proton ratio at BBN, changes from the stan- 
dard model value of this ratio will be reflected in its primordial abundance and, to a lesser extent, in the abundances of the 
other light elements produced during BBN. These relic abund ances therefore provide p robes of a universa l lepton asymme- 
tty. For further discussion and pr evious analys es, see e.g.. Wagoner. Fowler. & Hovle 1967; Reeves 1972; Yahil & Beauded 
1976; Beaudet & Gor el) 119761; iBeaudet & Yahil 19771; ISchramm & Steigman 19 79; Davi d & Reeves 1980; Scherrer 1983[ 
Freese et al.l 119831; iT erasawa & Sato 1985; Boesgaard & Steigman |1985l; iKan^ & Steigm anI 1 19921: iKohri. K awasaki & Sat3 
1997"; 'EsDosito et al. 2000; Ichikawa & Kawasaki 2003; Barger et al.| 120031: iKneUer & SteigmanI [20041: ICuoco et ah i2004t 
Serpico & Raffelt 2005; Steigman 2007; Popa & Vasile 2007, 2008. 

There is another way in which a significant neutrino degeneracy may influence the early evolution of the Universe. In the 
standard models of particle physics and cosmology, the energy density at, or just prior to, BBN is contributed by the cosmic 
background radiation photons, e^ pairs, and three flavors of extremely relativistic neutrinos. In the standard cosmology, the 
neutrino s constitute 40% of the to tal energy density. Any modification of the early-Universe energy density (or expansion rate; 
see e.g., ISimha & Steigmanll2008h can be parameterized in terms of the "effective" number of neutrinos by p ' = p+ ANt,/3^. 
For the standard models, N^ = 3 at BBN, while in the post annihilation epoch probed by the CMB, N^ = 3.046 (Mangan o et al.l 
5005) . A secondary effect of neutrino degeneracy is to increase the energy density in relativistic neutrinos predicted by the 
standard model, so that Nj^ N^ + I]qANi,(^q,), where 
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In general, if the three active neutrino flavors mix only with each other, neutrino oscillations ensure that their degeneracies 
equilibrate prior to BBN (Lunardini & Smirnov 2001; Dolgov et al. 2002; Wong 20 02;.Abazaiian. Beacom. & Bell.2002,) . In the 
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^ Charge conservation ensures that the very small proton-antiproton asymmetry, of order rjB , is balanced by a coirespondingly small asymmetry between 
electrons and positrons but, there are no such constraints on the size of any asymmetry between neutrinos and antineutrinos. 
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following we will assume this is the case and use = = ^t- As a result, 
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and 



Notice that, if = = ^r, then for < 0.2, SQAN,y(^Q) < 0.05. Later, in §6, we will relax this assumption so that 7''^^ 
and/or ^r)]- The only effect of non-zero values of and/or is to change the effective value of N^, while non-zero values of 
also contribute to N^, and, more importantly, they modify the neutron-to-proton ratio at BBN. 

2. NEUTRINO DEGENERACY AND BBN 

The stage is being set for BBN when the Universe is < 0. 1 seconds old and the temperature is > a few MeV, at which time the 
neutral-current weak interactions (e"*" + e" <-> i^ct + ; a = e, /i, r) are sufficiently rapid to maintain the neutrinos in equilibrium 
with the photon-e^ plasma. When the temperature drops below ~ 2 MeV, the neutrinos begin to decouple from the photon- 
plasma. However, the electron neutrinos continue to interact with the neutrons and protons through their charged-current weak 
interactions {n+Ve ^ p+e~ , p + v^ n+e^ ,n ^ p+e' + D^), enabling the neutron-to-proton ratio to track its equilibrium value of 
n/p= exp(- Am /kT), where Am is the neutron-proton mass difference. In the presence of a neutrino-antineutrino asymmetry for 
the electron neutrinos, where = £,(i'e) = -£,ii^e), the neutron-to-proton equilibrium ratio is shifted to n/p= exp{-Am/kT-£^e) 
which, for > results in fewer neutrons and less "^He. 

While the BBN-predicted abundances of D, ''He, and ^Li are most sensitive to the baryon density, that of '^He is very sensitive 
to the neutron abundance when BBN begins and, therefore, to any e-neutrino degeneracy. The primordial abundances of D, 
^He, or ^Li are baryometers, constraining t^b, while the ''He primordial mass fraction (Yp) is a leptometer, constraining . In the 
standard model of particle physics and cosmology with three species of neutrinos and their respective antineutrinos, the primordial 
element abundances depend on only one free parameter, the baryon density parameter, the post-e^ annihilation ratio (by number) 
of baryons to photons, 77B = n^/n^. This parameter is related to ba ryon mass de nsity parameter, Q,^, the present-Universe ratio 
of the baryon mass density to the critical mass-energy density (see lSteigmanll2006l) by 

7710= 10'%B/n^ = 273.9 r^B/i^ (5) 

The abundance of ''He is sensitive to (^e and to the early universe expansion rate S = H' /H = (p' /pY^^ = (1 + 7ANi//43)'''^ 
([Steigman 2007). In contrast, the abundance of ''He is relatively insensitive to the baryon d ensity since, to first order, a l l neutrons 
availa ble when BBN begins are rapidly converted to ''He. To a very good approximation (iKneller & Steigmanll2004t ISteigmaiil 
l2007h . 

Yp = (0.2482 ±0.0006) + 0.0016(r;He- 6), (6) 

where 

77He = mo+ 100(5- l)-5746/4. (7) 

In eg. [6l the effect of incomplete neutrino decoupling on the ^He mass fraction is accounted for according to the results 
of lMangano et al.1 (l2005l) . 

In contrast to ^He, since the primordial abundances of D, ''He and ^Li are set by the competition between two body production 
and destruction rates, they are more sensitive to the baryon density than to a non-zero lepton asymmetry or to a non-standard ex- 
pansion rate. For a primordial ratio of D to H by number, yu = 10^(D /H)p, in the range, 2 < ^ 4, to a very good approximation 
(IKneller & Steigmanll2004tlSteigmanll2007h . 

= 2.64(1 ±0.03)(6/77d)'', (8) 

where 

7/D = 7/10-6(5- l) + 5e./4. (9) 

The effect of incomplete neutrino decoupling on this prediction is at the ^ 0.3% level (iMangano et al.ll2b05h . about ten times 

smaller than the overall error estimate. 

A g ood fit to the primordial ratio of ^Li to H by number, = 10'"(Li/H)p is provided by (IKneller & Steigmanl2004tlSteigmaril 



l200^ . 



= 4.24(1 ±0.1)(77Li/6)^ (10) 

where 

7?u = ??io- 3(5-1) -76/4. (11) 

In our previous paper (ISimha & Steigmanll2008 h we assumed that S,e = and concentrated on the constraints on 7710 and 
N,y which follow from BBN and the CMB (supplemented by large scale structure (LSS) data). Here we will first set Nj^ = 3 
(S = 1) and repeat our analysis for 7710 and 6. For consistency, in this case we will need to confirm that |6| is bounded to 
be sufficiently small to justify the assumption that N^ = 3 + J^aA^i^(£,a) ~ 3. Then, we will relax the assumption that Nj^ = 3 
and use a combined BBN/CMB/LSS analysis to constrain all three parameters. In this case, any contribution to J^aA^ij{£,a) is 
automatically accounted for in our self-consistent determination of N,^. 

Before discussing our results, the primordial abundances adopted for our analysis are outlined. 
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Fig. 1 . — In the left hand panel the 68% and 95% contours in the § - '»yio plane from BBN and the adopted relic abundances of D and ^He. In the right hand 
panel, the BBN contraints convolved with the CMB/LSS constraint on 1710. 

3. ADOPTED PRIMORDIAL ABUNDANCES 

Since deuterium is destroyed as gas is cycled through stars, the deuterium abundance inferred from high redshift (i.e. young), 
low metallicity (i.e. little stellar processing) QSO absorption line systems should provide the best estima te of its primordial abun- 
dance. T he weighted mean of the seven, high red s hift, low metallicity D / H ratios fro m Buries & Tytler (1 1998ah . [Buries & Tvtleii 
il998b), Pettini & B owed (1200 Ih . ICMeara etaP (1200 ll) . iKirkmanetaP (l2003h . Wl/leara et alJ (120060 and, most recently, from 
[Pettini et al.. (.2008.) is 
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For this relic abundance. 
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Since the post-BBN evolution of He is also monotonic, with its mass fraction, Y, increasing along with increasing metallicity, 
at low metallicity, the "^He abundance should approach its primordial value Yp. Observations of helium and hydrogen recom- 
bination lines from low-metallicity, extragalactic Hii regions are of most value in determining Yp. At present, corrections for 
systematic unce rtainties do minate the esti mates of the observationally-inferred "^He primordial mass fraction and, especially, its 
error Following lSteigmanI (l2007l) : ISimha & Steigmanlf20(j8h . we adopt for our estimate here. 



Yp = 0.240 ±0.006. 



In this case. 



: 0.88 ±3.75. 



(14) 



(15) 



While the central value of Yp adopted here is low, the conservatively-estimated uncertainty is relatively large (some ten times 
larger than the uncertainty in the BBN-predicted abundance for a fixed baryon density). In this context, it should be noted that 
although very careful studies of the systematic errors in very limited samples of H ii regions provide poor estimators of Yp as a 
result of their uncertain and/or model-dependent ex t rapola t ion to zero metallicity, they ar e of value in providing a ro bust upper 
bound to Yp . Us i ng the results oflOhye & Skillmanl (l2004 . lFukugida & Kawasaki I (l2006h . and lPeimbert et all (l2007h . we follow 
ISteigmad ( l2007h : [Simha & SteigmanI ( |2008h in adopting a ^ 2cr upper bound of Yp <0.255. This upper bound to Yp corresponds 
to an upper bound to rjus < 10.25. 

Although the BBN-predicted relic abundance of lithium provides a potentially sensitive baryometer ((Li/H) oc ij^^, for rjn) > 4), 
its post-BBN evolution is complex and model-dependent. In addition, lithium is only observed in the Galaxy, in its oldest, 
most metal-poor stars in galactic globular clusters and in the halo. However, these oldest galactic stars have had the most time 
to modify their surface lithium abundances, leading to the possibility that the observationally-inferred abundances may require 
large, uncertain, model-dependent corrections in order to use them to infer the primordial abundance of ^Li. In the absence of 
corrections for depletion, dilution, or gravitational settling, the data of Rvan et al. ( 2000.) and .Asplund et al.. (2004) suggest 

[Li]p = 12 + log(Li/H) = 2.1±0.1. (16) 

In contrast, in an attempt to correct for evolution of the surface lithium abundances. lKorn et al.l (l2006l) use their observations of a 
small, selected sample of stars in the globular cluster NGC6397, along with stellar evolution models which include the effect of 
gravitational settling to infer 

[Li]p = 2.5±0.1. (17) 

In the following analysis the primordial abundances of D and "^He adopted here are used, along with CMB/LSS data, to 
estimate ryio, £,e, and AN^. Then, using the inferred best values and uncertainties in these three parameters, the corresponding 
BBN-predicted abundance of ^Li is derived and compared to its observationally inferred value. 

4. NEUTRINO DEGENERACY ^0); STANDARD EXPANSION RATE (N„ = 3) 

If attention is restricted to BBN alone, then there are two parameters inferred from observables, {f^o, ?7He}, and two unknowns, 
{vw,ie}- These are related by 

e.=4(77D-?7He)/579 (18) 
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Fig. 2. — In the left hand panel the 68% and 95% contours in the ^e— plane inferred from BBN and the adopted relic abundances of D and ''He, along with 
the CMB constraint on rjiQ. In the right hand panel, the BBN contraints convolved with the CMB constraints on and »yio. 

and, 

r/io = (574ryD + 577He)/579. (19) 

In the left hand panel of Figure [T] are shown the 68% and 95% contours in the - ijio plane derived from the adopted BBN 
abundances. Notice that and 7710 are virtually uncorrected. From this approach it is found that = 0.035 ± 0.026 and, that 
r/io = 5.88;!;y 33. The CMB/LSS provide an independent constraint on the baryon density, r/io = 6.14;|;J5 |[ (Simha & Steigman 2008). 
Since these two estimates for r^io are in agreement, we may use the CMB/LSS data to further restrict the allowed parameter space 
in the - 7710 plane, as shown in the right hand panel of Figure[T] In this case, the value of is unchanged and 7710 = 6.1 l^o [ [. 

Using these results for {r?io,Cc}' the primordial abundance of ^Li can be predicted using equations (10) and (11). In this case 
we find that rju = 6.05;!^ [j and [Li] = 2.63^o Q5. The non-zero value of which is consistent with BBN and the CMB is incapable 
of reconciling the BBN-predicted and observationally inferred primordial abundances of ^Li. 

An alternate approach which mixes BBN and the CMB provides another way to constrain the lepton asymmetry. Since the ^He 
abundance is much more sensitive to than is the deuterium abundance, equation (7) can be used along with the CMB value 
of 7710 to constrain ^g. This results in a virtually identical bound to that found from BBN alone, = 0.037±0.026. In contrast, 
combining the deuterium abundance with the CMB value of 7710 using equation (9) gives = -0.176±0.272, a range that is too 
broad to be of much value. 

From these results we see that, at 95% confidence, j^^ | < 0.09, so that I]QANiy(^Q) < 0.01. Such a small neutrino degeneracy 
has a negligible effect on the universal expansion rate during radiation domination, confirming the validity of the assumption that 
S = I (N,^ = 3). The corresponding lepton asymmetry t/l = 0.072 ±0.053, while "small", is orders of magnitude larger than the 
universal baryon asymmetry (t^b = 10~'*'77io = 6 x 10"'*'). 

5. NEUTRINO DEGENERACY ^0); NON-STANDARD EXPANSION RATE (N,. =/3) 

In this case, with N^ and ^j, free, there are three parameters to be determined: {7710, Nj^, ^g}, but only two useful relic abundances 
{yDjYp}. To constrain these three parameters it is necessary to employ the CMB/LSS data along with that from BBN. To obtain 
the best constraints a bit of care is required. While the CMB/LSS pr ovide a very good co nstraint on the baryon density parameter, 
they are less successful, at present, in constraining N^. For example. lSimha & SteigmanI (i2008.) find from the CMB/LSS data that 
while 7710 = 6.14+[5 [^, N^ = 2.9!^|| (or, S- 1 = -0.008!j][568)- Such a large rang e in will dilute t he constraint on ^g, resulting 
in a less than optimal constraint on given the available data. However, as ISimha & SteigmanI (12008) have noted, from the 
CMB/LSS the constraints on 7710 and N^ are virtually uncorrelated. Therefore, in the following we adopt for our observational 
input {77i(),77D,77He} and derive N,^ and from them. Since ^g is invisible to the CMB (except due to its contribution to the 
radiation energy density, which is accounted for in N^), we may combine our constraints on the {Niy,^g} combination with the 
independent constraint on N,^ from the CMB, to further constrain this combination. 

With S (Ny) and ^g to be constrained from 7/10, rjo, and rjug, equations (7) and (9) may be recast as 

S- 1 = (5797710 -5747/D-57yHe)/2944 (20) 

and, 

6 = (IO67710- 1007/D-67/He)/736. (21) 

In the left hand panel of Figure |2]are shown the 68% and 95% contours in the ^g - plane derived from the adopted BBN 
abundances in combination with the constraint on rjio from the CMB/LSS. From this approach we find ^g = 0.073;!;o Q57 and N^ = 
3.7±0.9. 



The CMB/LSS provide an independent constraint on Ni, , N,y = 2.9;I;o g (Simha & Steigman 2008). Since these two estimates 
are in agreement, the CMB/LSS data may be used to further restrict the allowed parameter space in the ^g - Nj^ plane, as shown 
in the right hand panel of Figure|2] In this case, we obtain ^g = 0.056±0.046 and N^, = 3.3;|;o g. 

Using these results for {7;io,Niy,^g}, the primordial abundance of ^Li can be predicted using equations (10) and (11). In this 
case we find [Li] = 2.62'^^-^. A non-zero value of fg along with a non-standard expansion rate, consistent with BBN and the 
CMB, are stiU incapable of reconciling the BBN-predicted and observationally inferred primordial abundances of ^Li. 
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Fig. 3. — The left hand panel shows the probability distribution for , inferred from BBN and the adopted relic abundances of D and ''He, for the case where 
N,y = 3. In the right hand panel, the dashed curve is the BBN constraint convolved with the CMB/LSS constraint on tjiq alone, and the solid curve is the BBN 
constraint convolved with the CMB/LSS constraint on t][o and N^ . 




Fig. 4. — The dashed curve shows the probability distribution for N^ for = inferred from BBN and the relic abundances of D and ''He convolved with the 
CMB/LSS constraints on rjio and Nj, . The solid curve shows the same for ^f. ^0. 

An alternate approach, mixing BBN and the CMB, provides another way to constrain the lepton asymmetry. Since the "^He 
abundance is much more sensitive to than is the deuterium abundance, equation (7) can be used along with the CMB/LSS 
values of 7710 and to constrain This results in a virtually identical bound, = 0.053±0.046. In contrast, combining the 
deuterium abundance with the CMB value of ryio using equation (9) yields = -0.061;!;q 37^, whose range that is too broad to be 
of much value. 

From these results, the effect of the neutrino degeneracy on the universal expansion rate during radiation domination may be 
computed, yielding Yia^^A^a) < 0.03 at 95% confidence. The corresponding lepton asymmetry, ryL = 0.115 ±0.095, while 



constrained to be "small", may nonetheless be orders of magnitude larger than the universal baryon asymmetry (t^b ■ 
6 X 10""'). 
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6. DISCUSSION AND CONCLUSIONS 

The left hand panel of Figure |3] shows the probability distribution of for the standard expansion rate S = (N^ = 3), 
derived after marginalizing over r^io , from BBN and the adopted primordial abundances of D and "^He. The right hand panel of 
Figure|3]shows the probability distribution of for the more general case where a non-standard expansion rate S 7^' (Nj^ ^ 3) is 
allowed. The constraints are based on combining BBN and the adopted primordial abundances of D and "^He with the independent 
constraints from the CMB/LSS. Additional information from the CMB/LSS constraints on r/io are used before marginalizing over 
r/io and N^ to produce the dotted curve, while additional information from the CMB/LSS constraints on both ryio and Nj^ are used 
before marginalizing over r^io and N,y to produce the solid curve. 

The dashed curve in Figure|4]shows the probability distribution of N^, for = 0, derived from BBN and the adopted primordial 
abundances of D and "^He with the independent constraints from the CMB/LSS after marginalizing over 7710 . The solid curve 
in Figure |4] shows the probability distribution of N,y for the more general case where a non-zero neutrino degeneracy 
is allowed. The constraints are based on combining BBN and the adopted primordial abundances of D and "^He with the 
independent constraints from the CMB/LSS. Additional information from the CMB/LSS constraints on 7710 and N^ are used 
before marginalizing over 7710 and ^g . Both of these constraints on are consistent with each other and with the standard model 
prediction of N^ = 3. 
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Our results can be used to constrain any deviation in the universal expansion rate from its standard model value due to the 
increase in radiation energy density from neutrino degeneracy. Our constraint on the neutrino degeneracy leads to Sq, AN,y(^Q.) < 
0.03 at 95% confidence. In addition, the constraint on the neutrino degeneracy yields a corresponding constraint on any lepton 
asymmetry, ?7l: 77l = 0.072 ± 0.053 for N^, = 3 and r/L = 0. 115 ± 0.095 when N^, 7^ 3. 

Using the constraints on 7710, N^,, and , the BBN-predicted primordial abundance of ^Li may be inferred. For N^y = 3, [Li] = 
2.63;I;0 og ^^'^ fo'" 7^3' [Li] = 2.62;!;!! n^- Both of these are considerably higher than the value ([Li]p = 2.1 ±0.1) determined from 
observations of metal-poor halo stars (iRvan et al.l(l2000h . lAsplund et al.l(l2006h) with out any correction for depletion, destruction, 
or gravitational settling. If, however, the correction proposed by iKorn et al.l (l2006h is applied, the predicted and observed ^Li 
abundances may, perhaps, be reconciled. It remains an open question whether this lithium problem is best resolved by a better 
understanding of stellar physics or, if it is providing a hint of new physics beyond the standard model. 

In our analysis we have assumed that = £,fi = S,t- Suppose instead that = £,t ¥ or, that = Cm t^C^ llDolgov et alj 
l2004l) . Since our constraints on come from BBN, and they constrain to be sufficiently small that the allowed degeneracy 
has minimal effect on the universal expansion rate (AN^(Cp) < 0.01), the only way to probe non-zero values of £^^=£,11= or 
£,a = £,T, is through their effect on the expansion rate (S or, AN^). In these cases, it is possible that £_a may be ^ £,e- For Nj^ = 
3.3;!^g, I]QANiy(Cct) < 1.7 at 2a. If it is assumed that = Cm ~ ^^^^ ^ 2. 34 and |?7l| ^ 5.0. If, instead, it is assumed 
that Cm = Ce ^ Cr (or, vice-versa for C^^), then |Cr | ^ 4.12 and |77l| ^ 7.6. 

Of course, our results a re sensitive to the relic abundances we have adopted. The simple but accurate fitting formulae 
dKneller & SteigmanlllOOl we have used make it easy to reevaluate our constraints for any adopted abundances. The con- 
straint on Cc is sensitive to the "*He abundance and is relatively insensitive to small changes in the D abundance. For example, 
we repeated our analysis for a different primordial '^He mass fraction, Yp = 0.247 ± 0.004. This alternate abundance, in com- 
bination with the D abundance used in this paper, and the independent constraints on 7710 and N^, from the CMB/LSS from 
ISimha & SteigmanI (|2008) yields £, = 0.023±0.041. 

Our results here are similar to, but considerably more restrictive than those of iBarger et al.l (l2003h and of ISerpico & RaffeE 
(|2b05), due to the improved constraints on Ni, and 7710 from t he WMAP 5 - year and oth er CMB and LSS data. The analysis 
described here seems to be related to that in recent papers by Popa & Vasild J20071 12008 ^. However, we fail to understand how 
they derive their constraints and why they find so much tighter bounds on Ni, and so much weaker bounds on Ce. 

Except from its contribution to the radiation energy density and the early Universe expansion rate, a lepton asymmetry in the 
neutrino sector is invisible to the CMB. Future CMB experiments will impr ove the constraint on by measuring the neutrino 
anisotropic stress more accurately. According to lBashinsky & SeljakI (l2004l) . PLANCK should determine N^, to an accuracy of 
(t{Ni,) ^ 0.24 and CMB POL, a satellite based polarization experiment, might improve it further to a{Niy) ^ 0.09, independent of 
the BBN constraints. Although still too large to provide a measure of the neutrino degeneracy, the tighter constraint on Nj^ can 
be used to further narrow the allowed ranges of Ce and N,y shown in Figure |2] 
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